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Abstract 

For any perturbative series that is known to fe-subleading orders of perturbation theory, we utilise the 
process-appropriate renormalization-group (RG) equation in order to obtain all-orders summation of series terms 
proportional to a n log n_fe (/x 2 ) with k = {0, 1, 2, 3}, corresponding to the summation to all orders of the leading 
and subsequent-three-subleading logarithmic contributions to the full perturbative series. These methods are 
applied to the perturbative series for semileptonic fa-decays in both MS and pole-mass schemes, and they result 
in RG-summed series for the decay rates which exhibit greatly reduced sensitivity to the renormalization scale 
fi. Such summation via RG-methods of all logarithms accessible from known series terms is also applied to 
perturbative QCD series for vector- and scalar-current correlation functions, the perturbative static potential 
function, the (single- doublet standard-model) Higgs decay amplitude into two gluons, as well as the Higgs- 
mediated high-energy cross-section for W + W~ — * ZZ scattering. The resulting RG-summed expressions are 
also found to be much less sensitive to the renormalization scale than the original series for these processes. 

1 Introduction 

The renormalization group equation (RGE) has long proven useful as a means of improving and extending results 
obtained from perturbative quantum field theory. In addition to giving rise to scale-dependent running parame- 
ters (coupling constants and masses) and concomitant scale properties (e.g. asymptotic freedom), the RGE can 
also be utilised to determine scale-dependent portions of higher-order contributions to perturbative expressions. 
For example, if the two-loop contribution to a physical process has been determined via explicit computation of 
pertinent Feynman diagrams, the RGE then determines all leading-log and next-to-leading contributions to all sub- 
sequent orders of perturbation theory. We denote such logarithms to be "RG-accessible." In the present paper we 
demonstrate how closed-form summation of such RG-accessible logarithm contributions is obtained for a number of 
physical processes whose field-theoretical series are known to two or more nonleading orders of perturbation theory. 
Consider a perturbative series of the form 

oo oo n 

S[xQi),LQi)] = J2 z" 5 " \ xL \ = J2^2 T n,kX n L k ( L1 ) 

n=0 n=0 k=0 

occurring within a physical decay rate T or measurable cross-section a, where x(n) is the running coupling constant 
(for QCD x(/j,) = a s (y)/ir) and where L(/j,) is a logarithm regulated by the renormalization mass scale y, that may 
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or may not also depend on a running mass: 

L^) = \og(^\ (1.2) 
If m is a running mass, then 

A* 2 ^r = rnj m [x(iJ,)] = -m[^x + -fix 2 + j 2 x 3 + ...] (1.3) 



If m is a pole mass (or for scattering processes, a kinematic variable), then j m as defined by (1.3) is zero. 

For example, in the fully MS expression for the semileptonic b — > u£~Di rate obtained from five active flavours, 
to is the running mass mb(p), 

r = |^ KMf5[^),£M] (i.4) 



and the successive-order series coefficients within S[x,L], as defined by (1.1), are |l| 

T ,o = 1, Ti,o = 4.25360, T X)1 = 5, T 2:0 = 26.7848, T 2)1 = 36.9902, T 2<2 = 17.2917 (1.5) 

The five active-flavour pole-mass expression for the same rate is obtained by replacing m,b(p) with the renormalization- 
scale independent pole mass mP h ole in (1.4) and (1.2), as well as a concomitant alteration of the following series 
coefficients jjj 

Ti l0 = -2.41307, T 1A = 0, T 2 , = -21.2955, T 2 ,i = -4.62505, T 2 , 2 = 0. (1.6) 

Suppose for a given scattering or decay process that the series S[x, L] is known to some order of perturbation 
theory: 

S NL = T ,o + (Ti, + T U L) a: (1.7) 

gNNL = S NL + fo Q + ^ + ^tfj x 2 (1 g) 

= gNNL + (Tg Q + T3 ii + Ta + T3 3l3) ^3 (1 9) 

5 iV 4 i = ^JV 3 L + ^ o + ^ iL + ^ + T4 sL 3 + TaaL ^ x \ ( L10 ) 

These next-to-leading (NL) and higher-order expressions exhibit scale dependence as the magnitude of L increases. 
However, higher order polynomial coefficients of L can be determined via an appropriate RGE. For example, in 
b — ► u£~i?£ the application of the RGE to the known two-loop (NNL) MS expression for the rate is sufficient to 
determine the three-loop coefficients T 3 , 3 , T 3j2 , and T 3A : for n f = 5, T 3j3 = 50.914, T 3<2 = 178.76, and T 3)1 = 249.59 
0]. This procedure is taken a step further in ref. 0, in which the four loop coefficients T^a, and T4 j2 are 

determined via the RGE for this same process. Estimates for T 3 arc also seen to determine T4.1, yielding an S N L 
expression characterized by only two unknown coefficients (T3 and I40) whose parameter space can be limited 
by the constraint that successive orders of perturbation theory decrease in magnitude 

In the present work, we wish to show how all RG-accessible logarithms may be summed if S is known to a given 
order. Specifically, we shall obtain explicit all-orders summations for the following four series, as defined by the 
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intermediate expression in ( 



1.1): 



5 [x(ii)L(ji)] = T ,o + T 1A xL + T 2 , 2 x 2 L 2 + T 3 ^x 3 L 3 + ... = £] T n;n x n L n (1.12) 

oo 

51 = Ti, + T 2A xL + T 3 , 2 x 2 L 2 + ... = £ T^x^L)"" 1 (1.13) 

n=l 

oo 

5 2 = T 2 , + Tg^xL + T^ 2 x 2 L 2 + ... = £ T„, rl _ 2 (xL)"- 2 (1.14) 

oo 

5 3 [a:(/i)L(^)] = T 3 , + T^xL + T 5 , 2 a: 2 i 2 + . . . = £ T nj „_ 3 (xL)"- 3 (1.15) 



n— 3 



The appropriate RGE ^/i 2 (r or a) = 0^ is seen to determine all series coefficients of S„ in terms of its leading 

coefficient T n ^, thereby facilitating the construction of RG-summed (RGE) perturbative expressions to any given 
order of perturbation theory: 

Srgtz = Sq[xL] + xSi[xL] (1.16) 
Srge = s o[xL] + xSi[xL] + x 2 S 2 [xL] (1.17) 
S Rck = S [xL] + xSi[xL] + x 2 S 2 [xL] + x 3 S 3 [xL] (1.18) 

These RGE expressions are seen to exhibit reduced sensitivity to the renormalization scale /i even when the 
logarithms L are quite large. Compared with the truncated perturbative series, these resummed expressions more 
effectively implement the underlying idea behind the RGE, namely that the exact (all-orders) expression for any 
physical quantity is necessarily independent of the scale-parameter \i. 

Although RGE determinations of higher-order terms have been known for some time to be of value in extracting 
divergent parts of bare parameters Q , the principle of incorporating all higher-order RG-accessible terms available 
to a given Feynman-diagram order of perturbation theory was, to the best of our knowledge, first articulated by 
Maxwell |^] as a method for eliminating unphysical renormalization-scale dependence. The all-orders summation of 
leading logarithms has been subsequently applied by Maxwell and Mirjalili || to moments of QCD leptoproduction 
structure functions and to NNL-order correlation functions. Such a summation of leading-logarithm contributions to 
all orders has also been explicitly performed by McKeon to extract one-loop RG functions from the effective actions 
of </> 4 -field theory in four dimensions and 6 -field theories in three dimensions . In Section ^ of the present work, 
we extend McKeon's summation procedure to derive closed-form expressions for all-orders summations of leading- 
( 1.12 ), NL-(1.12), NNL- (1.14), and N 3 L-logarithms (1.15) by using the RGE appropriate to the perturbative series 
flLll) within the QCD expression for the inclusive semileptonic B-decay rate. Such summations enable one to 
construct RGE perturbati ve e xpressions inclusive of up-to-three nonleading logarithmic contributions to all orders 



of the perturbative series (1.1). 

In Section ||, these results are applied to the b — > u£~p£ rate computed to NNL order by van Ritbergen 0, 
later extended via Pade-approximant methods to a subsequent ./V 3 L prediction 0. The renormalization-scale 
dependence of the unsummed perturbative rate truncated to a given order is shown to be much greater than that 
of the RGE rates obtained from the same perturbative expression. 

In Section |J, RGE expressions are obtained for the decays b — ► u£~D( and b — > cl~vi in the "pole mass" 
scheme in which only the couplant a s (ii)/n exhibits renormalization-scale dependence. This scheme, already 
known to have difficulties for the b — > u case [Q, exhibits a rate which increases with the renormalization scale fi, 
making the identification of a "correct" or optimal value of /i problematical. However, RG-summation is shown 
effectively to remove such /i-dependence, leading to reliable order-by-order pole-mass-scheme predictions for the 
b —t u semileptonic rate consistent with the b — * u rate obtained from an MS-scheme inclusive of a running b-quark 
mass. RGE-expressions are also obtained for the b — > c semileptonic rate based upon its (approximately-) known 
NNL series [pi and its Pade-estimated N 3 L series in the pole-mass scheme M. 
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The RGE appropriate for the perturbative series for semileptonic B-decays in the pole mass scheme is also the 
appropriate RGE for the fermionic vector-current correlation function utilised to obtain QCD corrections to the 
cross-section ratio er(e + e~ — > hadrons)/er(e + e~ — > In Section ^| we obtain RG-summation expressions for 

the QCD series embedded within the vector-current correlation function that include all higher-order logarithmic 
contributions that are accessible from the three fully- known nonleading orders of perturbative corrections in the MS 
scheme. We are thus able to compare directly the renormalization-scale dependence of the unsummed series S NL 



fll.Tp, S NNL M, and S N " L (M) to their corresponding RGS expressions S^s (MgD. S«ge (Mjlli and ^ags 



(1.18). We find that the latter expressions provide a set of virtually scale- independent order-by-order perturbative 
predictions for the vector correlator. 

In Section ^, we show how the use of process-appropriate RGE's can be used to obtain RGE perturbative 
expressions for a number of other processes. We obtain full RG-summations for 

1. the momentum-space series for perturbative contributions to the QCD static-potential function, 

2. the gluonic scalar-current correlation function characterising scalar gluonium states in QCD sum rules, 

3. the (standard-model-) Higgs-mediated cross-section W£W£ — > ZlZl at high energies, which is characterised 
by the two physical scale parameters s and Mb, 

4. the decay of a standard-model Higgs boson into two gluons [a process also characterised by two physical scales 
{Mh and M t ) in addition to the renormalization scale /i], and 

5. the fermionic scalar-current correlation function that characterises both Higgs — > bb decays and scalar- meson- 
channel QCD sum rules. 

We also discuss how RG-summation of the two scalar-current correlators considered in Section || removes much 
of the unphysical /^-dependence of the unsummed series at low s that would otherwise percolate through QCD 
sum-rule integrals sensitive to the low-s region. 

In Section 7 we summarise our paper. We discuss not only the reduction of /i-dependence via RG-summation, 
but also the comparison of RGS results with those of unsummed series when minimal sensitivity or fastest apparent 
convergence is used to extract an optimal value for the renormalization scale. 

Finally, an alternative all-orders summation procedure to that of Section || is presented in an appendix. 



2 RG- Summation of Logarithms for Semileptonic B-Decays 



For semileptonic b-decays, the //-sensitive portion of the rate ( |l.4| ) must, as a physically measurable quantity, 
exhibit renormalization scale-invariance: 

M 2 ^{K(^] 5 5[x( M ),L( M )]} =0. (2.1) 
This constraint is easily seen to lead to the RGE 

[l-2 7m (x)}^-+(3(x)^ + 5 7m S = 0, (2.2) 



where 



/3(x) = M 2 — x( M ) = -(f3 x 2 +pix 3 +l3 2 x 4 + ...) , (2.3) 
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[x(n) = a s ( /j)/ 7r] and where the anomalous mass dimension is the series defined by (1.3). Substitution of the series 
expansion (1.1) into the RGE yields the following series equation: 



= (l + 2 l0 x + 2 7 i.x 2 + 2 72 .x 3 + . . . ) J2 J2 Tn,kkx n L k ~ 1 

n=0 k=l 
oc n 

- ((3 x 2 + f3 lX 3 + f3 2 x 4 + ..•)££ T nik nx n - l L k 

oo n 

- 5 ( 7o x + 71a: 2 + 72a; 3 + ■ ■ ■ ) mi T ^x n L k . 



(2.4) 



n=0 k=0 



2.1 Evaluation of So 



To evaluate 5o[a;L], as defined by ( 1.12 ), we use (2.4) to extract the aggregate coefficient of x n L n 1 and to obtain 
the recursion formula (n > 1) 



nT n , n - [f3 (n - 1) + 57 ] T n -i, n -i = 0. 
We multiply (|2.5|) by and sum from rt = 1 to oo to obtain the differential equation, 



(1-^)^M-57o5oN = 0, 
du 



(2.5) 



(2.6) 



where Sq[u] is given by ( |l.l6| ) with xL replaced by u. The solution of ( |2.6|) for the initial condition So[0] = T is 

5oH=To,o(1-/3om)- 57o/ *. (2.7) 

For the special case of pole- mass renormalization schemes [y m [x] = 0], So = To o = 1, corresponding to the complete 
absence of x n L n terms from the perturbative series (1.1) when n > 1. 



2.2 Evaluation of Si 



To evaluate Si [u] , as defined by ( |OJ ) with u = x(p,)L(n) we first extract the aggregate coefficient of a;™L™ 2 from 
the RGE (|J) for n > 2: 



=(n - 1)T„, „_i + 27o(n - l)T„_i, „_i - f3 a {n - l)T n _ Xi „_ 2 
- 0l(n— 2)T„_ 2 , n _2 - 57 T!„_i ! „_ 2 - 57iT„_ 2 , n-2 ■ 



If one multiplies ( |2.8| ) by u n and then sums from n = 2 to infinity, one obtains the differential equation 

(1 - f3 u)^± - (/3 + 5 7o )SiM = 5 7l 5 [u] + (Pm - 2 7o )^. 

du ait 

We find it convenient to re-express this equation in terms of the variable 

w = 1 — (3qu 
and the constant 
57o 



A = 



Pa 



We see from $L7§ that if T , = 1, then 
So = w~ A 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 
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and find from (2.9) the following differential equation for Si: 
dSi 1 + A 



dw 



-S 1 = Bw 



-A-l 



Cw 



-A-2 



where 



B = (Afa - 5 7l ) /po 
C = A(2 l0 -p 1 /(3 ). 

For initial condition Si\ u —o = Si\ w —i = 2i ; o, the solution to ( [2.13 ) is 

Si = Bw~ A + [Ti >0 - B + C log(w)] w-^ 1 

with w, A, B and C respectively given by ( gig ), ( {TO] ), ( gig ) and (gl|). 



2.3 Evaluation of 5*2 

The aggregate coefficient of x n L n ~ 3 in cq. (|2.4|) is (n > 3) 



=(n - 2)T„,„_ 2 + 27o(n - 2)T„_i i „_ 2 + 271 (n - 2)T„_ 2 , „_2 - Po(n - l)T„_i, „_ 3 

- 0i(n - 2)T„_ 2 , n-3 - /?2(n - 3)T„_ 3! „_ 3 - 57 T n _i, „_ 3 - 57iT„_ 2 , n-3 - 57 2 T„_ 3! „_ 3 



If one multiplies ( |2.17j ) by u n 3 and sums from n = 3 to infinity one finds from the definitions 

00 

S [U] = l + J2 T n,nU n , 

00 

S 1 [u] = Y, T n,n-lU n - 1 , 
n=l 
00 



n=2 



[following from (p,.12| -L14)1 that 



dS 2 (2/? + 5 7 q) 5 = (Au-2to) dgi | Q3 2 u-2 7l ) dS | (/3i+5 7l ) g | 5 72 g 
du 1 — (3qu 2 1 — /9qu dtt 1 — /3qu du 1 — /Jqw 1 1 — Pqu 



(2.13) 



(2.14) 
(2.15) 



(2.16) 



If we incorporate the change-of-variable ( 2.10| ) in conjunction with the solutions ( 2. 12] ) and (2.16) for So 
respectively, we find that 



dS 2 , (2 + A) 



S2 = Dw 



-A-l 



Ew 



-A-2 



+ Fw~ A ~ 2 log(w) + Gw A 3 + Hw A 3 log(w) , 



-A-3 



.-A-3 - 



dw w 

where the constants {D, E, F, G, H} are given by 
D = [pxAB + f3 2 A - (fa + 5 7l ) B - 5 72 ] /#, 
£=(270-^) AB + [(1 + A) (Ti,o-S)-C]| 
F = [Afa ~ 5 7l ) C/0Q 
G=[(l + A) (T lfi -B)-C\ (270-^) 

F=(2 7 o-|)(l + A)C 



2 11 -^)A+(B- Ti, ) (ft + 5 7 i) //3 
Po, 



(2.17) 

(2.18) 
(2.19) 
(2.20) 

(2.21) 
and Si, 

(2.22) 

(2.23) 
(2.24) 
(2.25) 
(2.26) 

(2.27) 
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with constants {A,B,C} given by (2.11), (2.14) and (2.15). The solution to the differential equation (2.21) with 
initial condition ftL^o = ftlu^i = Tb.o is 



S 2 =-w 



+ Gw- A - 2 
2.4 Evaluation of S3 



A + {E- F)w~ A ~ 1 + Fw~ A ^ logH + (r 2t0 ^ F ) ,r 1 - 

H 



(2.28) 



log(iy) + -J w ~ lo S 



The aggregate coefficient of x n L n in (2.4) is 

=(n - 3) [T„,„_3 + 27 T„_i i „_ 3 + 2jiT n - 2 , n -s + 272^-3, n-3] 

- ft(n - l)T n _ Xi „_ 4 - ft(n - 2)T„_ 2 , «-4 - th(n- 3)T„_ 3 , „_ 4 - ft(n - 4)T„_ 4 , „_ 4 

— 5 7 oTn-i, „_ 4 — 57iT„_2, n -4 — 57 2 T„_ 3 ^ „_ 4 — 57 3 T , „_ 4j „_ 4 . 

To evaluate the series 

00 

S3M = X! T ".™- 3M 



(2.29) 



n— 3 



(2.30) 



n— 3 



WO multiply by u n 4 , sum from n = 4 to infinity, and, as before, make the eq. ( [2.10| ) change of variable 

w = 1 — ftit. We then find that 



dS 3 3 + A 

-: 1 03 —Kw 

aw w 



■Mw~ A ~ 2 + Nw~ A - 2] 



log(w) +Pw A 3 + QtiT A-3 log(w) 



A-3 1 



-A-3i„„2, 



iog^(w) + tz^r- 4 - 4 + ^ur A ~ 4 logH + y«j^- 4 \og\w) 



.-A-i 



-A-4i2, 



(2.31) 



by utilizing the explicit sol ution s fl2.12p , ( |2.16| ) and ( ]2.28|) for {5 , ft, ft}, as defined by (|2.18[) , (|2.19[ ) and ( |2.20| ). 
The new constants within ( 2.31| ) are 



K=^-(p 3 + Bfc + Dp 1 /2)-±- 
Po ft 



573 + (572 + ft) B + (5 7 i + 2ft) - 



M = 



9 ft 
272 

Po 



2 7l - ^ I B 1 2-,, - — ) - 



ft\ D 



ft 



,4 



ft 



+ [(57a + ft) (B ~ T ll0 ) + (5 7 i + 2ft) (F - E)] 

+ [(T 1 . Q -B)(1 + A)~C}^ + [E(1 + A)-F(2 + A)]^ 

ft ft 



AT = {(Aft - 5 72 ) C + [(A - l)ft - 5 7 i] F} 

ft 
ft, 
(57i + 2ft) 



ft 

P = ( 2- , - -f - ) ( 1 + . 1 ) ( T\ Xl -B) - C] + (270 - h ) [, l + . 1 ,£ - , 2 + . 1 )F 



ft 



T 2 , - - - S + P - ^ 



ft 



CV - (2 .1) j r 2 , - y -E + F 



(2.32) 
(2.33) 

(2.34) 
(2.35) 
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Q = 



u 
v 

Y 



2 7i " w ) C + ( 2 70 - ] F 



Pa 



(1 + A) - [(571 + 2ft) G + (P - (2 + A)G)ft] — 



if 



P = (/M - 5 7 l) 



(2 + A) ( T 2i0 - - - E + F ) - G 



27o-g ][(2 + A)G-fl] 



2 7 o 



(2.36) 
(2.37) 
(2.38) 
(2.39) 
(2.40) 



The constants {A, B, G, D, E, F, G, if} within ( f^3l| ) are respectively given by {(|]TT|), j2J^ ), ( pTl5| ) , ( |2~23| - |2~27| )}. 
The solution to ( [2.31 ), subject to the initial condition Sz\ w= i = P^o, is 



K 



5,3 =^w~ A 



-A-l 



log(to) + (P - Q + 2i?)w- A - 2 + (Q - 2P)uT A ^ log(iu) 



Rw 



-A-2 



log 2 H+(-f -f + f -P G 27? Pi.u ) 1 ' 



Pur A ' 3 log(ui) + — ur A_3 log 2 (u;) + — w~ A ~ 3 '~ ; 



log 3 (u) 



(2.41) 



where w = 1 — Pqu as in ( 2. 10| ) . 



3 Semileptonic b ^ u£ vi Decay in the MS Scheme 



In this section, we consider the semileptonic decay b — > u£~9g in the MS scheme. The decay rate is given by (1.4) 
in terms of the series (1.1). The coefficients in this series are fully known to two loop order and a re g iven by (1.5). 
The logarithms L(ijl) within the series are characterized by a running 6-quark mass, as given by (1.2) and (1.3). 



The coefficients (1.5) are listed for five active flavours, appropriate to analysis in an energy region containing 
TTib(mb). Consequently, the running &-quark mass rnf,(/z) and the running couplant x(fx) = a s (/z)/7r should be 
characterized by rif = 5 values for the RG functions 7 m [x] and f3[x\: 



7o 



253 

71 = "72~' 72 = 
23 29 

A) = pi = — , 02 



7.41986, 73 = 11.0343 
9769 



12' 



12' 



3456 



, #, = 18.8522 



(3.1) 
(3-2) 



Given the computed values of Po and Pjojjl | , it is straightforward to calculate the NL and NNL RG-summations 



for the series S, as defined in ( 1.16|) a nd (]l.l7 ). The constants {A, B,C, . . . , H} that characterise the summations 
So, Si and S2 are obtained via (|t|), ( |3.2|) and (1.5) from their definitions in Section Ej: 



A =—, B = -™**, <7=™°, D = 26.4461, 
23 3174' 529 

E = - 58.8224, F = - 3 -^, G = 25.5973, H = 1439220 



(3.3) 



279841 



279841 
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Eqs. (2.12), ( 2.16j ) and (2.28) then lead to the following closed- form expressions for the summations So, Si and S2' 

-60/23 



So 
Si 



23 
~ 12 

18655 



23 



3174 



1 



12 



-60/23 



f ^?log 



S 9 =13.2231 



529 



23 

1 ~ \2 X ^ L ^ 



23 
12 

-60/23 



x(fj,)L(fx) 



23 
12 



x(fi)L((x) 



-83/23 



(3.4) 



(3.5) 



3171350 
47 - 4897+ ^7984T l0g 

61.0515 + 25.5973 log 



23 
23 

1 ~ \2 X ^ L ^ 



1 ~ % X ^ L ^ 



-83/23 



719610 , 2 

log 

279841 5 



23 

1 ~ \2 X ^ L ^ 



23 

1 ~ i2 X ^ L ^ 



We first wish to compare the ^(-dependence of the 2-loop order expression 
yNNL 



K. 



[m fe (//)] 5 [1 + (4.25360 + 5L(//))x(/x) + (26.7848 + 36.9902L(^) + 17.2917L 2 (Ai)) x 2 {/j)] 



for the reduced rate (JC = Gj?\V u b\ 2 / 192ir 3 ) to that of the corresponding RG-summed expression 



-nNNL 
1 RGT, 

JC 



= K( M )] 5 [So + Six(ji) + S 2 x 2 (n)] 



-106/23 



(3.6) 



(3.7) 



(3.8) 



with So, S± and S2 given by (3.4), ( |3.5| ) and ( p.6[ ). To make this comparison, we evolve the running coupling and 
mass from initial values x(4.17GeV) = 0.0715492 and mb(4.17GeV) = 4.17 GeV ||, where the former value arises 
from rif = 5 evolution of the running coupling from an assumed anchoring value x (Mz) = 0.118000/7T |l(J, and 
where the latter value is the n/ = 5 central value in ref. []Tl"f for mb(nib). Thus x(p), to&(/x) and are fully 

determined via (1.2) and the RG-equations ( jl.3| ) and (|1|), with j m - and /3- function coefficients given by ( |3.l| ) and 
dp). 

In Figure 0, we use the u/ = 5 evolution of x(/i) and m&(/x), as described above, to compare r w,vi (3.7) to 
Tags' (3.S). It is clear from the figure that T^q-^ is almost perfectly fiat. By contrast, the naive rate r 7V N L is 
strikingly dependent on the renormalization scale fi, and does not exhibit any local extremum point of minimal 
sensitivity. Thus, RG-summation of leading, next-to-leading and next-to-next-to-leading logarithms is seen to 
remove the substantial theoretical uncertainty associated with the choice of /i from the (fully known) two-loop 
order b — > ut~T>g rate. 

It is useful to examine how the reduced b — > ut~vi rate develops in successive orders of perturbation theory. 
For example, the one-loop rates 



-— = [m b {n)f [1 + (4.25360 + 5L(ji))x(ji)] 

JC 

y<NL 

-^p = ko*)] 8 [So + s x x{n)] 



(3.9) 
(3.10) 



can be compared to the corresponding higher precision results of eq. (3/7) and (3J3). Three-loop order (N 3 L) 
reduced rates can be estimated through incorporation of an asymptotic Pade-approximant prediction of the three- 
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loop coefficient T3 = 206 The three-loop order expression for the reduced rate [] 



^N 3 L 



K 



= [m b (^i)f\ 1 + (4.25360 + 5L(/x))a;(/x) + (26.7848 + 36.9902L(/i) + 17.2917L 2 (/i)) x 2 ^) 



(206 + 249.592L(/z) + 178.755L 2 (//) + 50.9144L 3 (/i)) x 3 (fi) 



(3.11) 



can then be compared to its RG-summation version 



^N 3 L 
R.GT, 

K 



= [m b (/x)] 5 [S + Six(fi) + S 2 x 2 (fi) + S 3 x 3 (p)] 



(3.12) 



with So, Si and S2 respectively given by (|3.4|), (|3.5| ) and ( |3.6| ). The RG-summation S3, is obtained via (2.41). 
Given the estimate = 206, the known values (3.3) of {A, B, . . . , H} and values of {K, M, N,P, . . . , Y} defined 
via (p2|)-(|2~40|), 



K = -14.3686, M = 146.729, N = 50.9925, P = -317.085, Q = -148.520, 
R = -15.1138, U = 189.048, V = 83.3941, Y = 8.75961, 



(3.13) 



we find that 



4.7895 



23.^)60/23 



+ 



[60.617 + 25.496 log(l - ffxL)] 
(1 - ff xL) 83 / 23 



(1 

[-198.79- 118.29 log(l - '§xL) - 15.11410^(1- '§xL)\ 
(1 - ff xL) 106 / 23 

[348.96 + 189.05 log(l - §xL) + 41.697 log 2 (l - §xL) + 2.9199 log 3 (l - §xL)] 

23 



(3.14) 



(1 - ff rrL) 129 / 23 

The bold- face number 348.96 is the only coefficient in the above expression dependent upon the asymptotic 
Pade-approximant estimate for T%fi- We have included this estimate in order to demonstrate RG-summation 
incorporati ng a three-loop diagrammatic contribution to X^o; when such a calculation is performed, the factor 
348.96 in (|3.14|) should then be replaced by T 3:0 + 142.96. 

We consider the /i-dependence of three non-leading orders of perturbation theory first for the case in which 
logarithms are not summed to all orders. Figure || displays a comparison of the "unsummed" one-, two- and 



three-loop order reduced rates T NL /1C, T NNL /IC and T n3l /1C, respectively given by @, (§7j]) and |y]]). The 



//-dependence of all three orders is evident from the figure. Such //.-dependence can be used to extract NL and 
N 3 L values for T via the minimal-sensitivity criterion of ref. ]l~2]] . Curiously, T NL and T N L are both seen to have 
comparable minimal-sensitivity extrema (1801 GeV 5 and 2085 GeV 5 ) at values of \x much less than mbim^). r NNL 
exhibits some flattening between these extrema (~ 1900 GeV 5 ) over the same range of /t, but with a continued 
negative slope. Indeed, one can employ fastest apparent convergence |l3|] to choose /t for Y NNL such that \T NNL (jj,) — 
T NL (fi)\ is a minimum, and to choose /i for T n3l such that|r Af3i (//) - T NNL (n)\ = 0. As evident from Fig. § 
the former criterion leads to a value for /t (2.85 GeV) quite close to that value at which T NL (n) has an extremum 
(2.7 GeV), corresponding to T NNL /JC = 1888 GeV 5 . The latter criterion indicates that T N ' iL should be evaluated 
at the point where T NNL and T N L cross, a point noted previously Q to be virtually indistinguishable from the 
minimal-sensitivity extremum for T N L (/i). 

The point we wish to make here, however, is that all such values extracted for /i differ substantially from 
mi(ra(,), in which case progressively large powers of large- logarithms L(ji) = log[/i 2 /m 2 (/i)] enter the successive 



expressions (3.G), (3.7) and (3.11) for the NL, NNL and N 3 L rate T. Moreover, in the absence of minimal-sensitivity 



1 Only the non-logarithmic three-loop coefficient 206 is estimated; the remaining three logarithmic coefficients in (3.11) arc obtained 
via RG-methods in ref. fcj . 
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or fastest-apparent-convergence criteria for extracting /k, even the N 3 L rate exhibits a ±8% spread of values over 
the range m(,/2 < fj, £ 2mf,. 

RG-summation eliminates renormalization scale-dependence as a cause of theoretical uncertainty. In Figure 
we compare RG-summed versions of the reduced rate r^ GS (3.10), r^Y^ (3.8) and r^ G ^ ( |3.12 ). These three 
rates exhibit virtually no /i dependence whatsoever; rather, RG-summation is seen to lead to clear order-by- 
order predictions of the rate that are insensitive to [i. We see from Fig. || that r^ GS //C = 1646 ± 2GeV 5 , 

r ;?G£ /£ = 1816 ± 6 Gev5 . and r «Gs/^ = 1912 ± 4 Gev5 over the ( more or less ) physical range of \i considered 
in Fig. 0. Theoretical uncertainty in the calculated rate is now almost entirely attributable to truncation of the 
perturbation series to known contributions, an error which is seen to diminish as the order of known contributions 
increases. 

It is important to realize, however, that these scale-independent predictions necessarily coincide with the L(/i) = 
predictions of the unsummed rates (3.9), (3.7) and (3.11): T N L and r^ G ^ equilibrate when w = 1 i.e., when 
i(/i) = 0. Thus, one can argue that the summation we have performed here of all RG-accessible logarithms 
supports the prescription of identifying as "physical" those perturbative results in which /x-sensitive logarithms are 
set equal to zero. We must nevertheless recognize the possibility that the /i-sensitivity of the unsummed rates, when 
exploited by minimal-sensitivity or fastest-apparent-convergence criteria, is capable of leading to more accurate 
order-by-order estimates of the true rate than corresponding scale-independent RGE rates. In comparing Figs. |^ 
and H, it is noteworthy that the 1801 GeV 5 extremum of T NL /K, the unsummed one-loop rate, is quite close to 
r^Y^/ZC, the RG-summed two-loop rate. Similarly, the 1888 GeV 5 fastest-apparent-convergence value of T NNL /K, 
the unsummed two-loop rate, is close to r^ G ^. This train of argument would suggest that the 0(2085 GeV 5 ) 
extremum (or fastest-apparent-convergence value) of the unsummed rate T N L /JC may be a more accurate estimate 
of the true rate than r^ G ^. Such an argument, however, requires substantiation via explicit three- and four-loop 
order calculations, computations which are not yet available. 



4 Application to Semileptonic B decays in the Pole-Mass Scheme 



In the pole-mass renormalization scheme, the mass m appeari ng i n logarithms (1^2) is independent of the renormal- 
ization mass-scale /i. Thus the coefficients 7^., as defined in (.L2) are all zero. The constants {A, B, C, D, E, F, H, 
K, M, N, Q, R, Y}, as defined in Section are all zero as well. The nonzero constants are 

(4.1) 

r r, „ (4.2) 

r : -7 M1 — -2T-, AI (4.3) 

(4.4) 



rp Pi 

Po 






-Ti,o (fo 


P 2 l\ 

Po) 


1 

lh 


' T Pi 

Po 


- 2T 2 , 


Pi 
Po 


Po 







and the corresponding RGE summations appearing within S^ GS , Sj[q£, and S^q^ [(|1.16|), (1.17) and (1.18)] are 



So = 1 

S 1 [xL]=T 1>0 /{l-p xL) 

S 2 [xL] = T 2 , (l - P xL)- 2 + G(l - p xL)- 2 log(l - f3 Q xL) 



S 3 [xL] = P{1 - p xL)- 2 + 



V 



{T 3 ,o - P) + (71og(l - p a xL) + - \og z (l - (3 xL) 



(1 - PoxL)- 



(4.5) 
(4.6) 
(4.7) 

(4.8) 



In this section, we apply the above results toward the decay b — > u£ vg and b —* ci vg. The former rate is known 
fully to two-loop order in the pole mass scheme, though the result is argued to be of limited phcnomenological 
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utility Q. The latter rate has been estimated within fairly narrow errors to two-loop order as well 
been extended to a three- loop order estimated rate via asymptotic Pade-approximate methods . 



and has 



4.1 Pole Scheme Semileptonic b — > ut ui Decay 



ut V£ rate in the pole mass scheme is given by substitution of known values of the 



The two-loop order b 

coefficients {Ti j0 , 12,0, T 2 ,i, T 2>2 }, as listed in (1.6) into the series (|l.4|), with rrib(fi) replaced by 



pole 



as noted 



earlier. Although the reliability of the pole-mass scheme for this process is suspect because of the proximity of 
a renormalon pole JjJ], we have plotted this series for a range of x(fi) = a s ([j,)/ir and a choice for m? h ° le that 
will facilitate comparison with phenomenology already obtained from the corresponding MS process. We choose 
rif — 5 active flavours in order to explore the /i-dependence of the NNL rate in a region in which \x is considerably 
larger than mF h ole . Corresponding results for four active flavours are easily obtainable as well. The evolution of 
x{n) = a s (fi)/n for five active flavours ultimately devolves from a s (M z ) = 0.118 and leads to the same rif — 5 
benchmark value x(4.17GeV) = 0.071549 as noted in Section ||. Similarly we employ a value for mP h ole — 4.7659 
consistent to two-loop order with our use of the running mass value mb(mb) = 4.17 GeV, as obtained from the 



rif 



5 relation between 



pole 



and TOh(mb) of refs. O, [L 



ml = 4.17 GeV 



-x(4.17) + 9. 27793x 2 (4. 17) 



(4.9) 



We then see from Fig. that the //-sensitive portion of the known two-loop rate in the pole mass scheme, 



S NNL (fi) = 1 - 2.41307cc(/x) 



-21.2955 - 4.62505 log. 



x\a) 



(4.10) 



is indeed highly scale dependent. Specifically, we see that S nnl (/j,) increases monotonically with fi without ex- 
hibiting an extremum identifiable with a "physical" point of minimal sensitivity . 
In Figure H we have also plotted the RG-summed version of the 2-loop rate 



S 



NNL 
RGY, 



1 + x(fl)Sl 













x(ji) log < 




+ x 2 ( P )S 2 


x(ji) log < 









(4.11) 



with x(n) and rrP h ° e as obtained above. The summations S*i and S 2 are obtained via ( |4.6| ) and (4.7) using the 
rif = 5 pole-mass scheme values T^q = —2.41307, T 2y o = —21.2955 Q and rif = 5 QCD /3-function coefficients 
/3q = 23/12, Pi — 29/12, and (3 2 — 9769/3456. It is evident from the figure that renormalization scale- depen dence 
is considerably reduced by the summation of all orders of leading and next-to-leading logarithms in ( [3.11 ) . The 
increase of S^q^ with increasing \x is minimal compared to that of S nnl (/j,), the unsummed expression. 



In Figure [s], the comparison between 



po/e\5 qNNL 



pole ir, 

region rrv h j 2 



[i < 2m^ ole . The crossing point between the two curves necessarily occurs when L(/i) = 0, 



and 



po/e\5 qNNL 



) Srge i s exhibited over the physically relevant 



% 

pole 



corresponding to \l = mF h . Since S^f' is insensitive to fi, this crossover supports the expectation discussed 
in the previous section that the "physical" NNL rate is S NNL with fj, chosen to make all logarithms vanish. We 
would prefer, however, to is an almost scale-independent formulation of the NNL rate, thereby 

obviating any need to define a physically appropriate value of /i to compute a meaningful two-loop order result. 
We also note that the asymptotic large-/x result we obtain for the "reduced rate" 



r NNL 
1 RGH 

K 



(m 



pole\5 qNNL 



1829 GeV 5 



p,=Mw 



(4.12) 



is surprising close to the 1817 GeV 5 MS two-loop order ("unsummed" NNL) estimate obtained at /i = m^) = 
4.17 GeV, indicative of the utility of the pole-mass scheme when leading and next-to-leading logarithms are summed 



12 



to all orders. In the absence of such summation the pole mass expression S NNL (fi) spans values for the reduced 
rate between 1420 GeV 5 and 2060 GeV 5 as /i increases from 1 GeV to Mw, reflecting the problems with the pole- 
mass scheme already noted in ref. By contrast, the RG-summed reduced rate varies only from 1774 GeV 5 to 
1829 GeV 5 over the same region of fx. 



4.2 Pole Scheme Semileptonic b — > c£ vt Decay 

The semileptonic decay of B into a charmed hadronic state is given by the following decay rate in the pole-mass 
renormalization scheme S : 



r(6 -> ci~v t ) = 



(4.13) 



In ( 4.13| ), mi, an d rn c are /^-invariant pole masses, L(/i) = log (/i 2 /mfcTO c ) , x(/i) = a s (p)/iT, and F(r) is the form 
factor 

F(r) = 1 - 8r - 12r 2 log(r) + 8r 3 - r 4 . (4.14) 
All sensitivity to the renormalization scale [i resides in the series S[x, L], which may be expressed in the usual form 
S[x, L] = l + (Ti )0 + T hl L) x + (T 2 , + T 2A L + T 2)2 L 2 ) x 2 + (T 3 , + r 3 ,iL + T 3 , 2 L 2 



T 3 , 3 L 3 ) x 3 



(4.15) 



For four active flavours, the perturbatively calculated coefficients of ( 4.15 ) are T\$ = —1.67 and the (partially- 
estimated) coefficient T 2i o = — 8.9(±0.3) fi 15|. Except for T^q, the remaining coefficients in ( 4.15[) are accessible 
from the RG equation 







d_ 

dL 



f3(x) 



d_ 

Ox 



S[x, L] 



(4.16) 



These coefficients are T„,„ = [n > 1], T 2 ,i = T h0 /3 a = -3.479, T 3 ,i = 2T 2 , /3o + T 1>0 Pi = -42.4(±1.3), and 
T 3-2 = Ti,o/3 2 = — 7.25 ||. An asymptotic Pade-approximant estimate of T 3i0 = — 50.1(±2.6) has also been 
obtained in ref. |J. Consequently, one may list three orders for the /i-dependent portion of the b — » ci~Dt rate: 



S NL [x(ji),L(ji)]= l-1.67z( M ), 
S nnl [x(ji),L(ji)] = 1 - 1.67x(p) 



-8.9- 3.479L(/x)]a; 2 (At), 



S N L [x(p),L(n)] = S NNL [x(n),L(n)} + [-50.1 - 42.4L(/*) - 7.25L 2 (/x)] x 3 (/x), 



(4.17) 
(4.18) 
(4.19) 



with the caveat that NNL and N 3 L expressions have increasing theoretical uncertainty arising from the (small) 
estimated error in T 2 o and concomitant error in the estimation of T 3 - 

As before, we will compare the \i dependence of ( 4.17j) , ( l.lgj ) and ( 4.19| ) to that of the corresponding RG- 
summed expressions 



S 



%&: = l + x( M )S 1 [x(n)L(fx)] 

iNNL 



SrSv = 1 + ^)Si [x(ji)L(jj.)] + x'(p)S 2 [x(fj,)L(fi)] 



Q N 3 L 



s 



NNL 



x 3 {^)S 3 [x(n)L^)], 



(4.20) 
(4.21) 
(4.22) 



in order to illustrate how RG-summation of higher logarithms affects the order-by-order renormalization-scale 
depend ence of a per turba tive series. Figu re ^ displays a plot of the ^t-sensitive portions of the decay rate considered 
to NL (|4.17|) , NNL (|4.18|) and N 3 L (|4.19|) orders. As in |J § , m c is assumed to be m b /3. We have chosen the pole 
mass mb to be 4.9 GeV consistent with phenomenological estimates [@J^| The couplant x(fi) = a s (ii)/n is chosen 

2 Such an estimate is slightly larger than that based upon ( |i.s| ), as ( |j;,s| ) is truncated after two-loop order. 
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to devolve from a s (m T ) — 0.33 |17j via four active flavours, where po = 25/12, /3± = 77/24, [3 2 = 21943/3456 and 
/?3 = 31.38745. These choices permit careful attention to the 1.5 GeV < [i < low-scale region anticipated to 
correspond to the physical rate, although we have chosen to extend the range of /i to ~ Im^ in Fig. |6| 

Figure |^ demonstrates that the rate expressions appear to progressively flatten with the inclusion of higher order 
corrections, but that the residual scale dependence of each order remains comparable to the difference between 
successive orders. Figure ^ also displays rates proportional to the corresponding RG-summed expressions (4.20), 
4.21| ) and ( 4.22| ) based upon the same phenomenological inputs. The expressions for Si, S 2 and S3 are given by 



(4.7) and (4.8). It is evident from Figure || that the scale dependence of RG-summed expressions to a given 



order is dramatically reduced from the scale dependence of the corresponding non-summed expressions. 



Thus S 



NL 

i?G£' ^RGY, 



s 



NNL 



and Sj^2 are effective scale-independent formulations of the one-, two- and three- loop 
perturbative series within the b — > c£~9g rate; once again the summation of progressively less-than-leading loga- 
rithms in the perturbative series is seen to remove the choice of renormalization-scale /1 as a source of theoretical 
uncertainty to any given order of perturbation theory. 



5 The Vector-Current Correlation Function and R(s) 



The imaginary part of the MS vector-current correlation function for massless quarks can be extracted from the 
Adler function p8[. This procedure is explicitly given in |19] and leads to an expression in the following form: 



1 



Im n„(s) 



J2Q}S [s(A0,log(/i a /*)]> 

/ 



where a;(/x) = a s (/x)/7r and s is the kinematic variable p 2 [i.e. the square of the invariant mass in e H 



The series S[x,L] appearing in (5.1) is fully-known to N 3 L order: 

S N ° L [x, L] = 1 + x + (T 2 , + T 2A L)x 2 + (T 3 , + T 3>1 L + T^ 2 L 2 ) 



(5.1) 
hadrons] . 

(5.2) 



The full series is, of course, identifiable with the generic series ( |l . l[ ) [or (4.15)] provided T .o — 2~i,o = 1 and 
T n , n = [n > 1]. Values for the remaining vector correlation function constants in (5.2) are tabulated in Table [j] 
for three, four and five flavours. The coefficients T 2 fl and T3 are obtained from the results of ref. [Ell, and are 
well known from the standard expression for perturbative contributions to R(s). The remaining coefficients 



T 9 



Po, T3 1 — 2poT2 + PI, ^3,2 



Pi 



(5.3) 



are easily determined from the renormalization scale-invariance of the vector-current correlation function (5.1) 



/x^S^m), \og{n 2 /s)] = 



S\x,L] = 



This equation, of course, can be interpreted to reflect the imperviousness of the physical quantity 



R(s) 



hadrons) 3 

■+ jU + // ) 47T 



(5.4) 



(5.5) 



cr(e+e 

to changes in the choice of QCD renormalization-scale /1 [ ^o| . 

Equation (^J) is just the RGE ( |2.2[ ) with 7 m (x) set equal to zero — precisely the same RG-equation as applicable 
to the pole-mass scheme semileptonic b decay rates considered in the previous section. Consequently, the RG- 
sum mat io n of th e sc ries S\ x, L ] within ( |o.l| ) involves the same series summations So, Si, S 2 and S3 as those given 
by (15), (16), (17) and (18). The (nonzero) constants G, P, U and V appearing in these equations are found 
in terms of /3-function coefficients /3 , Pi, §2 and series coefficients T 10 (= 1), T 2 , and T 3fi via eqs. (4.1) -(4.4) 



These constants are all tabulated in Table pj, and are seen to fully determine the 0(N 3 L) RG-summed version of 
the series S[x,L], 



Sr G j:[x, L} = 1 + xSi[xL] + x 2 S 2 [xL] + x 3 S 3 [xL}. 



(5.6) 



14 



where x = x(fi) and L = log (/i 2 /s). For example, if n/ = 5, we see from @), @ and and the Table | 

entries for G, P, U, V, T2,o and Tk.o that 



Sl[xL]= (1-S*L) (5J) 

r n 1.49024 -1.26087 log (l-ffxZ) / 

5 2 [xL] = — ^ i 2 -^ (5.8) 

(1 - f§ xL) 



, r .._ 0.115003 [-12.9196 - 5.14353 log (1 - §xL) + 1.58979 log 2 (l-§xL) j 

63 xL — — -j- H — -3 (5.9) 

[1-§xL) 



In Figures ^ and §1 we compare the /i-dependence of the unsummed (^^) and summed ( |5.6| ) expressions for 
5[x(/x), log(/z 2 /s)], with the choice s = (15GeV) 2 . The running coupling constant x(fj,) is assumed to evolve via 
the n/ = 5 (four-loop-order) /3-function from an initial value x(M z ) — 0.11800/-7r|] Although Figure fj] does show 
a flattening of the unsummed expressions upon incorporation of successively higher orders of perturbation theory 
[S NL , S NNL , S N L ], Figure || demonstrates that the corresponding RG-summed expressions are order-by-order 
much less dependent on the renormalization scale /i. In particular, the full N 3 L summed expression ( |5.6| ) exhibits 
virtually no dependence on fi, but is seen to maintain a constant value S^ G ^ = 1.05372 ± 0.00004 over the entire 
i/s/2 < (j, < 2yfs range of renormalization scale considered. By contrast the unsummed expression S N L o f (|5.2| ) is 
seen to increase (modestly) over this same range from 1.0525 to 1.0540. The point marked FAC in Figure |7| is the 
intersection of the unsummed expressions for S NNL and S N L . This point is the particular choice of fj, at which the 
0(x 3 ) contribution to (5.2) vanishes, i.e., the point of fastest apparent convergence (FAC). It is noteworthy that 



this FAC value for S" 3i is quite close to the RG-summation value S^ G ^, a result anticipated by Maxwell Q. In 
other words, if one were to use the FAC criterion to reduce the theoretical uncertainty arising from /j-dependence 
in the original expression (O) for S N L {n) over the \fs/2 < /i < 2y / s range considered, the specific value one 
would obtain [Sp AG — S N L (26 GeV) = 1.05402] corresponds very nearly to the RG-summation value extracted 



from (154 [S$& = 1-05372 ± 0.00004] 



This behaviour is not peculiar to the choice of s in Figs. [7] and || In Figures ^ and 10, we consider S N L , S N L 
and Sx G x for y/s = 30GeV and y/s — 45 GeV. The FAC point each figure occurs at the value of /i for which 
S NNL (fi) — S N l (/j), as discussed above. In both figures it is evident that S N L has substantially more variation 
with /i than Sr G ^, which is virtually independent of /i. Nevertheless, both figures also show that the FAC point of 
S N * L is very close to the Sf^-level value. For ^5 = 30 GeV [Fig. § sf A L G = S n3l (h = 52 GeV) = 1.04700 
and S% G % = 1.04689 ± 0.00002. For ,/s = 45 GeV [Fig. Sg A % = S N3L (fi = 77 GeV) = 1.04369, and 
s rg& = 1-04360 ±0.00001. 

Similarly, the PMS points in all three figures, corresponding to maxima of S N L (/i), are very near the FAC 
points and also quite close to the Sf^s leveL In Fi S- @ [V 5 = 15 GeV], S^ s = (23 GeV) = 1.05403; in Fig. 
I [vs = 30 GeV], sf A [ s = S N * L (46 GeV) = 1.04701; and in Fig. [l(] [yfi = 45 GeV], S$f s = S n3l (70 GeV) = 
1.04370. Of course, equality between S N L (/J,) and S^ G ^. is necessarily exact when L = 0, i.e., when /i 2 is chosen 



equal to s. This is indeed the prescription employed in the standard 1 1 prescription relating S to R(s) |l 



R(s) =l + x(Vi) +T 2 , x 2 (V^) +T 3:0 x 3 (^) , (5.10) 

where T2 i o and T3.0 are given in Table [l] for nj — {3,4,5}. The point here, however, is that this prescription is 
justified not by the /j-invariance of S n3l [x(/j), log (/j 2 /s)] , the truncated perturbative series, but by that of S^ G ^, 
the perturbative series incorporating the closed-form summation of all RG-accessible logarithms within higher- 
order terms. Moreover, the nearly /j-independent result S^ G ^ appears to be quite close to the result one would 
obtain from S N L [x(/x),log (a* 2 /s)] either by imposing FAC or PMS criteria to establish an optimal value of /j, 

3 For purposes of comparing the /^-dependence of S n3l and S^q^, we are assuming (as in Sections^ and^) there to be no uncertainty 
in the value of ct s (/i). 
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as graphically evident from Figs. |^ and [l^. That optimal value for /i, however, is not /i = y/s, but a substantially 
larger value of \i for each case considered. 

These results are not peculiar to the choice rif = 5. Using the Table [I] entries for the parameters T 20 , T30, G, P, U 
and V appearing in (4.7) and (4.8), we find for rif — 3 that 



Srg^\ x > L] — 1 + 



[1.63982 -f log (l-fxL)] 



„3 ( 3397 \ 



and for n/ = 4 that 



1-fzL 

a; 3 [-8.97333- 8.99096 log 


(l-f.L) 2 
(1 - f .xi) + log 2 (1 




(5.11) 


(1 


- hL? 






a; x 2 [1.52453 - 


ilog(l-f|xL)] +a ; 3 


( 121687 \ 
l 180000/ 




.t 3 [-11.0096- 7.06715 log 


(i-i-^) + iiiog 2 ( 


1 - f 2 xL)] 


(5.12) 



(1 - f 2 xL) 



In Figures [ll] and we display the expression ( |5.12| ) for S^ G ^ together with (unsummed nj = 4 expressions 
for) S N L and S NNL for y^s = 4GeV and y/s = 8GeV. To generate these figures, we evolve x(fi) using an initial 
condition a;(4.17GeV) = 0.0716218 appropriate for nf = 4 and obtained from the threshold matching conditions 
HI to the n/ = 5 running couplant x(4.17GeV) = 0.0715492 evolved from x{M z ) = 0.118/rr Juj. As in Figs. | 

are both close to the very-nearly constant S^ G ^ level, although the PMS and 
FAC values for /i are substantially larger than y/s in each case. 

To conclude, the primary result of interest is that closed-form summation of all RG-accessible logarithms to 
any given order of perturbation theory leads to expressions (e.g. Fig. ^) that order-by-order are substantially less 
scale-dependent than the corresponding truncated series (Fig. Q). The scale-independence of S^ G ^ supports the 
prescription of choosing [i = yfs in the unsummed series S n3l \ji 2 , log(/i 2 /s)], since the summed and unsummed 
series coincide at this value of [a. However, the unsummed series S n!>l still exhibits noticeable scale dependence. 
The use of FAC and PMS criteria to find an optimal value for fj, for S N L leads to values for this unsummed series 
which are quite close to its RG-summation value. 

6 Other Perturbative Applications 

6.1 Momentum Space QCD Static Potential 

The momentum space expression for the perturbative portion of the QCD static potential function 



d3 9 iff.? ( 16tt 2 



v(r)= 1 Wf eiqr[ -— ]w 



log 



,2 



q z 



(6.1) 



is given by the integrand series 

W[x, L] = x + (T 2 , + T 2 ,iL) x 2 + (T 3 , + T 3 ,iL + T 3 , 2 L 2 ) x 3 + . . . (6.2) 



where x = a s (/x)/7T, L = log (/i 2 /(f 2 ), and where the series coefficients within (6.2) are 22 

T 2 . = 31/12 - 5n//18, T 2A =(3 = 11/4- n//6, 

J-//324, T 3A — 24 1 / — Man// 1 i t uibf/ o-±, — mo 



T 3:0 = 28.5468 - 4.14714n/ + 25n 2 /324, r 3)1 = 247/12 - 229n//72 + 5n 2 /54, T 3j2 = 



(6.3) 



The function W[x, L] is shown in |23j to satisfy the same RG equation (5^4) as the semileptonic b — > u and 6 — ► c 



decay rate in the pole mass scheme. Consequently, the closed form summation of all RG-accessible logs is given 
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by eqs. (J4.6|)— (4.8 ) , with the constants G,P,U and V as given by eqs. (4.1)— (4.4). Since Tq.q, as defined by the 



generic series form (1.1), is zero for the series ( |6.2| ), the series Sq = J2^=o T n ,n{xL) n is trivially zero. Thus, the 



RG-summation of W[x, L] is 

W RG z = xSt[xL] + x 2 S 2 [xL] + x 3 S 3 [xL]. (6.4) 
Noting that T lfi = 1 in (O), we find for n f = {3, 4, 5} that 



W, 



f [i-jlog(l-§^)]+^ 3 [im] 



(6.5) 



cc 3 



W ' 



[9.62511 - log (1 - §xL) + f§ log 2 (1 - jxL)] 
(1 - §xLf 

'< , - 2 [i-Iiog(i-i^)]-- 3 [iii] 



,3 [13 . 8688 _ 77693 lQg (1 _ K xL) + 5929 lQg 2 (l _ g^j ( 6 ' 6 ) 

25, 



(1 - J^XL 



W. 
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Iog(l-|xL)]-a» [ii] 



(6.7) 



[18.1104 - f> log (1 - |s£) + g log 2 (1 - IxL)] 
(1-%xL) 3 



In Fig. [13] we plot both the truncated series W N NL \x, L], consisting of the terms explicitly appearing in (6.2), 
as well as the co rres ponding RG-summed series ( |6.5[) for the nt — 5 case with |g| = 4GeV. As before, we determine 
x(fi) from eq. (U) with the initial condition a s (4.17GeV)/7r = 0.071549 devolving from a s (M z ) = 0.118000. 
Although the truncated series varies only modestly [0.0856 > W NNL > 0.0795] over the range |g|/2 < // < 2|g|, 
the RG-summed series is seen to exhibit less than 10% of the unsummed series' variation over this same range of 
H [0.08295 > W BGT , > 0.08233]. 

6.2 Gluonic Scalar Correlation Function 

The imaginary part of the correlator for gluonic scalar currents 

My) = ^r^G^(y)G^(y) (6.8) 
enters QCD sum rules pertinent to scalar glueball properties p4| , and is given by (s = p 2 , x(fx) — a s (fi)/ir) 

Imn G ( S ) = Im |i J d 4 y e** < 0|Tj G (y)j G (0)|0 >| = [xQi), log ( M 2 / S )] (6.9) 

where 

00 n 

S[x, L] = 1 + J2J2 T n,mX n L m . (6.10) 

The leading coefficients T n ^ m within (|6.10|) can be extracted from a three-loop calculation by Chetyrkin, Kniehl 
and Steinhauser [^5[ [l9] and are tabulated in Table |2[ 
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The correlator ImllG(s) is RG-invariant: /j, 2 dImIlG(s)/dfi 2 = 0. Consequently the series ( 6.10| ) can be shown 
to satisfy the RG- equation 



d_ 

dL 



d | 2g(s) 
dx x 



S\x,L] = 0. 



(6.11) 



Upon substituting ( 6.1 0| ) into ( 6.1 1| ) , it is straightforward to show that the aggregate coefficients of x n L n 1 , x n U' 
and x n L n ~ 3 respectively vanish provided 



nT n , n - flo(n + l)T„_i, „_i = 0, 

(n - l)T„,„_i - P (n + 1)T„_ 1] „_ 2 ~ PinT n _ 2 , n-2 = 0, 

(n - 2)T„ i „_ 2 - /9o(n + l)T„_i, n-3 ~ /3inT„_ 2 , n-3 ~ /M™ _ !) r n-3 : n-3 = 0. 



(6.12) 
(6.13) 
(6.14) 



We e mplo y the definitions ( |2.18| ), (J2JJ) and ( |2.20[ ) for So, Si and S 2 . By multiplying ( J6.12D by 



Q6.13D by m 



n-2 



and (|sj1 by u"~ 3 and summing from n = 1, 2 and 3 respectively we obtain the following three linear differential 
equations 



(i - M 



dSo 
du 



(1 - P u)^ - 3f3 Q S 1 



(1 - fl>«) 



dS 2 
du 



2/3 ^o = 
4PoS 2 = Pi 



t dSo 
du 
dSj_ 
du 



2S 
■3Si 



u— h 26 

du 



Given the u = initial conditions So = 1, Si — Ti.o, S2 
RG-summed version of the series ( 6.10 ) to NNL order: 



(6.15) 
(6.16) 

(6.17) 

T 2 .o and setting u = xL, we obtain the following 



Srg^ =S [xL] + xS x [xL] + x 2 S 2 [xL] 



T h0 -^log(l-p Q xL) 



+ x* 



111 



- ■ 

00 



(1 - p xL) 2 
(t 2 ,o - 



(1 - foxL) 3 



(6.18) 



+ x 



20i 



2fe 
/3o 



+ if) log(1 " /?oxi) + ^(l - PoxL) 



(1 - p xL) 4 



where x = a s (/u)/7r and L — log(^ 2 /s), as before, and where the f3k are as defined in (|2.3[). In Figure |l4| we compare 
the /z-dependence of x 2 {^)S^q^ to that of the corresponding truncated series 



NNL r 



, L] = x 2 [1 + (Tlq + Tx.iL) : 



(T 2 ,o + T 2 iL + T 2 , 2 L 2 ) x 2 } 



(6.19) 



for the n/ = 3 case with s/s = 2GeV. The evolution of x(fi) is assumed to follow an rif = 3 /3-function with 
the initial condition a s (m T )/7r = 0.33/7T ]l7|j . As evident from the figure, the seve re /i-dependence of x 2 S NNL is 
considerably diminished by RG-summation. The RG-summed expression [eq. ( |6 .18 ) multiplied by £ 2 (a0] falls from 
0.056 to 0.041 as n increases from 1 GeV to 4GeV. By contrast, the unsummed expression ( 6.19 ) falls precipitously 
from 0.259 to 0.017, a factor of 15, over the same range of /i. This unphysical dependence on renormalization 
scale suggests that S NNL [as in ( 6. 19] )] be replaced by ( |6.18 ) within sum rule approaches to the lowest- lying scalar 
gluonium state. 



6.3 Cross-Section a(W£W£ -> Z L Z L ) 

The scattering of two longitudinal Ws into two longitudinal Z's is mediated by the Higgs particle of standard- 
model electroweak physics. Assuming a single Higgs particle (devolving from the single doublet responsible for 
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electroweak symmetry breaking), one finds the cross-section for this process at very high energies (s >> Mjy) to 
be 



a(s) 



8^ 
9s 



5 2 ( M )%(M),log(M 2 /s)L 



(6.20) 



where g(fi) — 6A^jg(/x)/167r 2 , the quartic scalar couplant of the single-doublet standard model,^ and where the 
series S is pfll 



OO 71 



S[g,L] = l + J2J2 T ^9 n L r ' 



(6.21) 



with g = g(/x),and L = log(/i 2 /s). The constants T n ^ m are fully known to NNL order [^6| 

Ho = -10.0, T M = -4, T 2:0 = 93.553 + | log (s/M^ ) , T 2jl = 68.667, T 2 , 2 = 12. 
The RG-invariance of the physical cross-section a implies that the series S satisfies the RG equation 



d_ 

dL 



g , W{g) 
dg g 



S[g,L}=0 



(6.22) 



(6.23) 



This is the same RG-equation as ( |6.11 ) characterizing the gluonic scalar correlator, but with the (3- function appro- 
priate for the single-Higgs-doublet MS quartic scalar couplant g(fi) p7|: 



fa 



(3{g) = -p g 2 - Pig 3 - fcg 4 ■ ■ ■ 
-2, 0t = y, (3 2 = -27.803 



(6.24) 



Consequently the RG-summed version of S is given by ( |6.18 ): 

1 g[-W+f log(l + 2$L)] -18.4141 5 2 



Srgs = 



(1 + 2 5 L)2 

2 



(1 + 2 5 L)3 



g 



111.967+ | log (^) 



iff log(l + 2 5 L) 



log 2 (l + 2 5 L) 



(6.25) 



(1 + 2gLf 



6.4 Higgs Decay H 
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Higher-order expressions for the decay of a Higgs boson into two gluons have been obtained and studied both 
outside |2£| and within |2^, |l| the context of a standard-model single-doublet Higgs field. In the limit Mfj < < 
4M 2 , Mb — the latter decay rate is of the form 



V2G f M 3 h /72tt\ x 2 (h)S Km), log (>7™?(a0) , log (M 2 H /M 2 )] 



(6.26) 



Capitalized masses (Mn,Mt) denote RG-invariant pole masses, whereas mt(p) is the running t-quark mass. The 
series S within ( 6.26| ) is of the generic form (1.1) with L = log (^ 2 /m 2 (/i)) , but the coefficients T n>m are now 
dependent upon the RG-invariant logarithm T = log (M^/M 2 ). Using six active flavours to accommodate the 
running of mt(/i), one can extract from ref. p9| the following two subleading orders of series coefficients within S 



To.o = 1, T hQ (T) 



12 6 ' 2 



T 2 (T) = 146.8912 - + ^T 2 , T 2 l( T) = ^ - ^T, T 2 2 = H 

2,ov ; 48 48 ,u ; 16 8 ' ' 16 



(6.27) 



4 Aj^jg-(/i) is perturbatively related to its on-mass-shell value G as discussed in |26 
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RG-invariance of the physical decay rate (/j, 2 dT /d/j, 2 = 0) implies the following RG-equation for the series S[x, L, T] 
within fl6.26| ): 



„ , „ 3 nl s d 20{x) 



S[x,L,T] = 0. 



The rif = 6 values for the MS ft and 7 m -funetions are 

0{x) = - (fox 2 + (3 lX 3 + ftx 4 + ...), ft = 7/4, ft = 13/8, ft = -65/128 
7m (a) = - (702: + 7l x 2 + 72a; 3 + ...), 7o = 1, 7i = 27/8, 72 = 4.83866. 



(6.28) 



(6.29) 



Upon substituting the series ft as described above, into ( 6.28| ), one finds that the net coefficients of x n L n 1 , x n L n 2 
and x n L n ~ z on the left-hand side of (6.28) vanish provided the following recursion relations are respectively upheld: 



nT n>n - ft(n + l)T n _i, n _i = 

(n - l)r n>n _i + 27 (n - l)T n _ 1)n _x - o (n + l)T n _ ljT ,_ 2 - 0\nT n _ 2 , n -2 = 
=(n - 2)T n> „_ 2 + 27o(n - 2)T„_i ! „_ 2 + 271(71 - 2)T„_ 2 ,n-2 
- ft(n + l)T n _i >n _3 - /3inT„_ 2 ,n-3 - ft(n - l)T n -3,n-3 ■ 



(6.30) 
(6.31) 

(6.32) 



By multiplying (3.30) by u n , (6.31) by u n 2 and (6.32) by u n 3 , and by summing each equation from n = 1,2 
and 3, respectively, we obtain the following linear differential equations for the summations So (2.18), ft (2.19) 
and ft (pfj|): 



(1-^)^0 _ 20o S = 
dit 



(1 - ftu)-^ - 3ftft = ft u— ^ + 2ft - 2 7o — - 



dti 

dft 



dS 



du 



du 
dft 



dS 



du 



du 
dft 



(1 - o u)—^ - 4ftft = Pi + 3ft - 2 70 -^ + ft it— ^ + 2S - 2 7 i— ^ 



dit 



dft 



dit 



dS Q 



du 



(6.33) 
(6.34) 

(6.35) 



Given the u = initial conditions So = 1, ft = T\fi(T), ft = 12,0 (T), one can solve for ft[ii], ft[u] and ft[u]. As 
before, all-orders summation of the RG - acce ssible logarithms within the series S[x, L, T] is no w pos sible, given the 
explicit form of T\.o{T) and T 2 ,o(T) in ( 6.27 ) and the explicit ft and 7-function coefficients in ( |3.2£ ). We thus find 
that 



-iNNL 
~>RGT, 



—Sq[xL] + xSi[xL] + x 2 S2[xL] 



1 , xj^-fT+f lo g (l-lxL)]+^x 2 



(6.36) 



x 2 [134.801 -^T+^T 



529 rp 2 , / 21675 345 



392 



^T)log(l-|,L) + filog 2 (l-I,L)] 



(1 - ixL) 



6.5 Fermionic Scalar Correlation Function 

The imaginary part of the RG-invariant correlator for the fermionic scalar current 
j s (y) = m*(y )¥(t/) 

is 



Im II(s) = Im 



d 4 y e 1 ^ < 0|Tj fl (y)j s (0)|0 > 



3s 



(6.37) 



(6.38) 
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where the series S[x, L] is of the form ( [L~l| ) and has been fully calculated to N 3 L order |32). For nj — {3, 4, 5} the 
series coefficients T„ m are tabulated for (n, m) < 3 in Table ^. This correlation function is relevant both for QCD 
sum-rule analyses of scalar mesons, a topic of past and present interest p3| , and for the decay of a single-doublet 
standard-model Higgs boson into a 66 pair |30). RG-invariance of the correlator (/i, 2 dlmll(s)/d/i 2 = 0) implies 
the following RG-equation for the series S[x,L] within (6.38) p9[: 



S\x,L] 



(6.39) 



where L = log (fj, 2 /s), x — a s (fi)/ir, j3(x) is the /3-function series (2J3) and J m ( x ) is the 7-function series (1.3). 
Substitution of ([Ll]) into ( 6.39| ) leads to the following recursion formulae for the elimination of terms proportional 
to x n L n ~ l , x n L n ~ 2 , x n L n ~ 3 , and x n L n - 4 : 



= nT ntn - p (n - l)T„_i iri _i - 2j T n -i,n-i 

= (n - l)T„ jn _i - /3o(n - l)T„_i >n _2 - 0i(n - 2)T n _2,„-2 - 27oT n _i i „_2 - 27iT„_ 2 ,n-2 
=(n - 2)T„, n _ 2 - /9o(n - l)T n _i jT ,_ 3 - f3x(n - 2)T„_ 2 ,«-3 - fii(n - 3)T„_ 3! „_ 3 

~ 270-^71-1,71—3 — 27iT n _2 >n _3 — 272T n _ 3 n _3 
=(n - 3)T n ,„_ 3 - /3 (n - l)T„_i,„_4 - /3i(n - 2)T , n _2,, l - 4 - /3 2 (n - 3)T , n _ 3i „_ 4 

— Pz{n — 4)T„_4 >n _4 — 27oT n _i iT[ _4 — 27iT n _ 2 ,n-4 — 27 2 T„_ 3:Il _4 — 27 3 T , „_4 ! „_ 4 

We follow the usual procedure of 

1. multiplying ( 6.40| ) by and summing from n = 1 to 00, 



(6.40) 
(6.41) 

(6.42) 
(6.43) 



2. multiplying (6.41) by u n ~ 2 and summing from n — 2 to 00, 

3. multiplying ( 6.42| ) by u™ -3 and summing from n — 3 to 00, and 

4. multiplying ( 6.43| ) by «" -4 and summing from n — 4 to 00. 

Using the definitions (|I|), fl2~2fj|) and ( pp| ) for {5 , Si, S 2 , S3}, we then obtain the following four linear 

differential equations for these summations: 



(l_ / 3 oU )_£_2 7o S = 0, S [0] = 1 

(1 - Po^-r 1 - (A) + 2 7 o)Si = ftu-j^ + 2 7l 5 , 5i[0] = T h0 
au au 

(1 - - (2/3 + 2 7o )S 2 = 0iu-^- + (3 2 u-^- + (0i + 27OS1 + 272S0, S 2 [0] = T 2 , 



(1 - p u)^ - (3/? + 2 7o )S 3 =/?i«^ + 

du dit du 



/3 3 u^ + (2/3i + 2 7l )S 2 
clu 



+ (A + 272)^1 + 2 73 S , S 3 [0] = T 3 , 



The solutions to these equations are 
S [u] = (l-f3 u)- A 
S 1 [u] = C 1 {l-f3 u)- A + 



a , [T lfi - C! + C 2 log(l - J3 u)} 



(1 - f3 u) 



A+l 



(6.44) 
(6.45) 
(6.46) 

(6.47) 

(6.48) 
(6.49) 
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S 2 [u] 



(1 - pou)- 



[D 2 - D 3 + D 3 log(l - p Q u)] 



(1 - Pou) 



A+l 



[T 2 ,o -^-D 2 + D 3 + Dj log(l - Pou) + ^ log 2 (l ~ Pqu)} 
(1 - [3 u) A + 2 

&M= T (i-A)«) + (i - M a + i 

, F 4 - F 5 + 2F 6 + (F 5 - 2F 6 ) log(l - ftu) + F 6 log 2 (l - o u)} 



(6.50) 



(i - M A + 2 

+ f-F i + F 5 -2F e + F 7 log(l - A)^) + f log 2 (l - A)") + f log 3 (l - P u)] 



(1 - p u) A +3 



where 



A = 2 10 Po 

C\ = 2Pao/Po - 27i//?o 
C 2 = -2/3 l70 //? 2 

D 1 = [Ci (A/3! - p x - 2 7l ) + - 2 72 ] //3 

£> 2 = [Ci (2 7i - 2Aft) - C 2 Pi - + Ti j0 - 2 7l )] 

D 3 = C 2 (A(3\ - 2 7l ) /p 

Di = pi [d -(A + l) (Ti, - Ci)] /A, 

D 5 = -C 2 /3i( J 4 + l)//3 

Fx = [Ci (Ap 2 - 2 72 - /%) + D x (Apx/2 - 7l - A) + A/3 3 - 2 73 ] /#, 



Po 



D 3 (2 7l - APi) + D 2 {Afc -Pi- 2 7l ) 



Dt 



AP X - C 2 p 2 



+C\ (2 72 - 2Ap 2 ) + Ti, (A/3 2 - 2 72 ) - A/3 3 



F 3 = [£> 3 - Pi ~ 27i) + C 2 (Aft - 2 72 )] //3 
1 



(T 2 , - £>i/2 -D 2 + D 3 ) (Aft - 2 7l ) - D 4 ft + £ 3 (A + 2)ft - £> 2 (A + l)ft 



+C 2 p 2 + d(A + l)p 2 - T lfi {A + l)p 2 



F 5 = [D 4 (Aft - 2 7l ) - £> 3 (A + l)ft - £> 5 ft -C 2 (A + l)p 2 ] /ft 
F 6 = D 5 (Aft/2- 7l )/ft 

F 7 = pt [D A -{A + 2) (T 2 , - Dx/2 -D 2 + D 3 )} /p 
F 8 =ft [D 5 -D 4 (A + 2)} /p 
F 9 = -ft£> 5 (A + 2)/(2ft,) 



(6.51) 

(6.52) 
(6.53) 
(6.54) 
(6.55) 
(6.56) 
(6.57) 
(6.58) 
(6.59) 
(6.60) 

(6.61) 

(6.62) 

(6.63) 

(6.64) 
(6.65) 
(6.66) 
(6.67) 
(6.68) 



The 0(N 3 L) RG-summation of the series S[x(/i), \og(fi 2 /s)] appearing in the correlator ( |S.38 ) is then found 
to be 



S 



N 3 L 
AGS 



=S 



log — 



x 3 (/z)5 3 



"(M)logF- 



(6.69) 



where the RG-summations So, Si, S 2 and S3 are given by (|6l8|)-(|6~5ll) with u = x(n) log(/i 2 /s). 
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In Figure n~q we compare the /x-dependence of the RG-summed scalar fermionic-current correlator, 



Imn 



RGT, 



m 2 {jj)S 



RGT, 



to that of the correlator (6.38) when truncated after its fully-known 0[x ) contributions, 

3 n 



lmn N L ~ m 2 



log ( E- 



(6.70) 



(6.71) 



over the range y/s/2 < fi < 2*Js with y/s — 2 GeV. The coefficients T„ jm appearing in (6.71) are given in Table 
H We choose to work in the \/s — 2 GeV rif — 3 regime appropriate for QCD sum rule applications, where the 
couplant x(n) is large. The evolution of x(fi) is assumed to proceed via the n/ = 3 four-loop /3-function with initial 
condition x{m T ) = 0.33/7T [ jl7| . The running mass m(/x) is normalized to 1 GeV at /i = m T to facilitate comparison 
of ( |6~70| ) to (| ,7l[ ). In Figure [l5|, the unsummed correlator is seen to achieve a sharp maximum near 1.5 GeV, 
followed by a precipitous fall as [i approaches 1 GeV from above. By contrast, the RG-summed correlator exhibits 
a much flatter profile, falling from 1.88 GeV 2 to 1.75 GeV 2 as /x increases from 1 to 4 GeV. As in Section 6.2 for 
the gluonic scalar-current correlator, these results indicate that even V r3 L-order expressions for the perturbative 
contribution to QCD correlation functions exhibit substantial /^-dependence. Such dependence, which we have 
shown to be largely eliminated via the RG-summation process, would otherwise percolate through QCD sum-rule 
integrals as spurious sensitivity on the theory side to the Borel parameter. The incorporation of RG-summcd 
correlators within QCD sum-rule extractions of lowest-lying scalar resonances is currently under investigation. 



7 Summary 

In this paper we have explicitly summed all RG-accessible logarithms within a number of perturbative processes 
known to at least two nonleading orders, a procedure originally advocated by Maxwell Q. As anticipated, we have 
found the dependence on the renormalization scale /x in every case examined to be considerably diminished over 
that of the original scries' known terms. 

In semileptonic b — > u decays in the fully MS scheme (Section ||), we observe the intriguing possibility that 
PMS /FAC criteria for the unsummed series truncated to a given order may anticipate the RG-summed series for the 
next order of perturbation theoryj^] This behaviour, however, is not evident in the other processes we consider. For 
the vector correlation function (Section |J), PMS/FAC criteria for the unsummed series truncated to a given order 
coincide closely with the RG-summed series for that same order, but do not anticipate the level of the next-order 
RG-summation. In the pole-mass scheme version of semi-leptonic B-decays (Section ^), PMS and FAC criteria 
do not appear applicable to the unsummed series, which monotonically increase with the renormalization scale 
\x. Indeed, one of the virtues of RG-summation is the sensible scale-independent results it provides for inclusive 
semileptonic B-decays in the pole mass scheme, a scheme whose unsummed expressions for b — > u are already 
known to be problematical Q . 

The /x-independence of RG-summation, particularly for the vector-current correlation function (Section |5|), is 
seen to justify the prescription of zeroing all logarithms by setting the renormalization scale to /i 2 equal to the 
kinematic variable s. This prescription necessarily equates the unsummed and RG-summed series, and, since the 
RG-summed series is virtually independent of /x, the zeroing of logarithms in the unsummed series equilibrates it 
to the flat RG-summation level we obtain. 

In Section |^, RG-summation is also applied to the perturbative contributions to the momentum-space QCD 
static potential, the decay-rate of a standard-model Higgs to two gluons, the Higgs-mediated cross-section a{WW — > 
ZZ), and to two scalar-current correlation functions. Examination of these last two quantities in the low-s region 
appropriate for QCD sum-rules suggests the utility of RG-summation is reducing the unphysical scale-dependence 
of the perturbative QCD contributions to the field-theoretical side of sum-rules in these channels. 

5 Although this result incorporates a Pade estimate for an RG-inaccessible 0(x' s ) coefficient, this estimate occurs in both the 
unsummed and RG-summed expression. 
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Appendix: An Alternative Closed-Form Summation Procedure 

The body of our paper has addressed the evaluation of S n [xL] = X^fclo T n +k,k (xL) , where the full perturbative 
series S(x, L) — J2n°=o ^ n \ x ^\ ^ however, also possible to group the terms within S(x, L), as defined by eq. 
(1.1), such that the dependence of each series term on x and L fully factorises: 



S(x, L) = ^R n {x)L r ' 

oo 

R n {x) = ^T k , n x k ■ 



(A.l) 
(A.2) 



k—n 



Let u s supp ose , for example, that the series S{x,L) satisfies the RGE (2.5) appropriate to 6-decays. By substituting 
( |A~l"l ) into (|2~2"1), we find that 



R n +i{x) 



1 



n+1 



1 



2 7 (x) - 1 



0(x)— +5 7 (a;) 
da; 



R n {x) 



(A.3) 



where we have relabeled the anomalous mass dimension 7 m (x) — > 7(2;) to avoid any misinterpretation of the label 
m as a subscript. If 



R n ( x )=exp(-j^01dx')p n (x) , 



(A.4) 



then the recursion relation (A.3) implies that 



Pn+l(x) 



1 



0(x) 



-Pn{x) 



n+1 ^27(2;) — 1 J da; 
If one defines x to be implicitly a function of y via the equation 

jL X (y) = gfr) 

dy 2 7 (x) — 1 ' 



then ( |A.5| ) simplifies to 

P n+1 (x(y)) = J—^- Pn ( x (y)) 



n + 1 dy 



(A.5) 



(A.6) 



(A.7) 



in which case 



n=0 



n=0 



1 „ d" 



(A. 
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This last result implies via (A. 4) that the series S(x, L) is fully determined by knowledge of the log-free summation 
Rq, i.e. that 



S(x, L) = exp 



x(y+L) 



/?(*') 



R (x(y + L)) 



where x{y) is defined implicitly by the constraint 



y = 



27(^0- i d , 
p>') 



obtained by integrating (A. 6). 



Using lowest-order expressions (3{x) = —(3qX and 7(2;) = —702;, we find from (A. 10) that 
' -i + 27 log(a;) J +K . 



Po V x 



If we set 



2 7o W^ 



K=^ log (2 7o ) 
Po 



we find from (A.ll) that 



We = exp 



270 



(A.9) 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



Eq. (A.13) is the defining relationship for the Lambert VK-function W [exp (— /3oJ//27o)], as discussed in ref. 
Since 



1 



2 l0 W [exp (-/3 o2 //2 7o )] 
in the approximation /3(x) = —(3qx 2 , 7 (x) = —70^, we then find from (|A.S| ) that 



Ty [exp (-/W2 7 o)] 
W[exp(-/3o(y + £)/2 7 o)] 



570 //9 



1 



2 7o ^[exp(-/3 ( 2 / + L)/27 )] 



(A.14) 



(A.15) 



where 



i?o(x) =^T M x fc =5(x,0) 



fc=0 



and where 



y = - -T- + ^ log (2702;) 

Poz Po 



(A.16) 



(A.17) 



For the RGE (4.16), corresponding to ( |2.2| ) with 7 m (x) chosen to be zero, the solution ( JA. S| ) still applies provided 
x(y) is defined implicitly via the constraint 



y = 



da' 

P>0 



(A.18) 



25 



In the approximation /3(x) = — (3qx 2 , one can choose 



z(v) = -77- 

Poy 



(A.19) 



in which case 



S(x, L) = R 



1 — (3qLx 



(A.20) 



where the function Rq is as defined via (A.I 
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rif = 3 


rif = 4 


71/ = 5 




9/4 


25/12 


23/12 


Pi 


4 


77/24 


29/12 


Pi 


3863/384 


21943/3456 


9769/3456 


T2.0 


1.63982 


1.52453 


1.40924 


22,1 


9/4 


25/12 


23/12 


2*3,0 


-10.2839 


-11.6856 


-12.8046 


2b, 1 


11.3792 


9.56054 


7.81875 


T3.2 


81/16 


625/144 


529/144 


G 


-16/9 


-77/50 


-29/23 


P 


3397 

2592 

-8.99096 


121687 
180000 

-7.06715 


17521 
152352 

-5.14353 


U 


V 


512/81 


5929/1250 


1682/529 



Table 1: Constants for determining the 0(N 3 L) RG-summation of the 7V 3 L-order vector-current correlation func- 
tion. 





71/ = 2 


71/ = 3 


nf = 4 


71/ = 5 


71/ = 6 


2i,o 


6919 
348 


659 
36 


4999 

300 


4123 
276 


367 
28 


2i,i 


29 
6 


9 
2 


25 
6 


23 
6 


7 

2 


2^2,0 


246.434 


197.515 


150.210 


104.499 


60.3685 


22,1 


7379 

48 


2105 
16 


1769 
16 


4355 

48 


1153 
16 


32,2 


841 

48 


243 
16 


625 
48 


529 
48 


147 
16 



Table 2: NNL-order series coefficients within the gluonic scalar current correlation function. 





71/ = 3 


n/=4 


n f = 5 


2b, 


1 


1 


1 


2i,o 


17/3 


17/3 


17/3 


2i,i 


2 


2 


2 


32,0 


31.8640 


30.5054 


29.1467 


22,1 


95/3 


274/9 


263/9 


22.2 


17/4 


49/12 


47/12 


23,0 


89.1564 


65.1980 


41.7576 


23,1 


297.596 


267.589 


238.381 


23,2 


229/2 


22547/216 


10225/108 


23,3 


221/24 


1813/216 


1645/216 


ft 


47.2280 


31.3874 


18.8522 


73 


44.2628 


27.3028 


11.0343 



Table 3: iV 3 L-order series coefficients within the fermionic scalar current correlation function, as calculated in fl32| |. 
Also listed are the four-loop /3-function j34j] and 7-function [ |35"| coefficients 03 and 73 required for the evaluation 
of the series S3 . 
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Figure 1: Comparison of the ncxt-to-ncxt-to- leading- (NNL-) order unsummed (dotted line) and RG-summed 
(solid line) decay rates T/JC for b — > u£~P( in the fully-MS scheme with five active flavours (n/ = 5). The quantity 
JC = G 2 F \V ub \ 2 /1927T 3 . 
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Figure 2: Comparison of unsummed b — > u£ v t decay rates in the fully-MS scheme (rif — 5) truncated after 
NL-order (solid line), NNL-order (dotted line), and N 3 L-order (dashed line). 
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Figure 3: Comparison of RG-summation expressions for the fully-MS b — > ul v t decay rate (n/ 
the NL (solid line), NNL (dotted line) and N 3 L (dashed line) perturbative series. 
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Figure 4: Comparison of the large-/x behaviour of the NNL unsummed (dotted line) and RG-summed (solid line) 
decay rates for b — > ul~vi within the RG-invariant pole-mass scheme with five active flavours. 
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Figure 5: Comparison of the NNL unsummed (dotted line) and RG-summcd (solid line) decay rates for b — > ui 



(rif = 5) within the pole-mass scheme over the range ml° le /2 < jj, < 2m^ 
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Figure 6: Comparison of unsummed and RG-summed b — ► c£ vi decay rates T/K (rif — 4) in the pole-mass 
scheme, where JC = G% \V cb \ 2 /192tt 3 . The curves representing the unsummed rates are labelled by NL, NNL and 
N 3 L indicating the order at which they are truncated. Similiarly, the RG-summed curves are labelled by SNL, 
SN 2 L and SN 3 L. 
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Figure 7: Comparison of the following n/ = 5 vector-current correlation- function series when y/s = 15 GeV: 
S NL (solid line), the MS perturbative series S [x (fi), log (^ 2 /s)] truncated after NL-order contributions; S NNL 

(dotted line), the same series truncated after NNL-order contributions; and S n3l (dashed line), the same series 
truncated after N 3 L-order contributions. At the intersection of S NNL with S N L , the N 3 L-order contribution to 
the perturbative series is zero, corresponding to the point of fastest apparent convergence (FAC). 
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Figure 8: Comparison of the following RG-summations of the perturbative MS series within the n/ = 5 vector- 
current correlation-function when ^fs = 15GcV: S^q^ (solid line), the summation based upon the perturbative 
series S [#(//), log (/z 2 /s)] truncated after NL-order contributions; S^q^ 1 (dotted line), the summation based upon 

the perturbative series truncated after NNL-ordcr contributions; and Srq^ (dashed line), the summation based 
upon the perturbative series truncated after N 3 L-order contributions. 
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Figure 9: Comparison of S NNL (dotted line), S N L (dashed line), and S^ G ^ (solid line) expressions for the series 
S [#(//), log (^i 2 /s)] within the n/ = 5 vector-current correlation-function when ^fs = 30 GeV. 
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Figure 10: Comparison of S NNL (dotted line), S N L (dashed line), and S^ G ^ (solid line) expressions for the series 
S [x(fj), log (^i 2 /s)] within the n/ = 5 vector-current correlation- function when ^fs = 45 GeV. 
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Figure 1 1 : Comparison of S NNL (dotted line) , S N L (dashed line) , and S^ G ^ (solid line) expressions for the series 
S [x(/x),log (li 2 /s)] within the n/ = 4 vector-current correlation-function when = 4 GeV. 
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Figure 12: Comparison of S NNL (dotted line), S N L (dashed line), and S^ G ^ (solid line) expressions for the series 
S [x(fj,) 7 log (li 2 /s)] within the n/ = 4 vector-current correlation-function when = 8 GeV. 
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Figure 13: The momentum-space n/ = 5 static-potential function series W NNL (dotted line), and the correspond- 
ing RG-summation W^ L (solid line) with \q\ = 4 GeV. 
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Figure 14: The imaginary part of the gluonic scalar-current correlation function ( |6.9| ), as obtained from x 2 (n) S 
(dotted line) and from x 2 (/j,) (solid line) with nf = 3 and = 2 GeV. 
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Figure 15: The imaginary part of the fermionic scalar-current correlation function ( 6.38| ), as obtained from 
m 2 (^) S NNL (dotted line) and from m 2 {p) S%££ (solid line) with n f = 3 and ^ = 2 GeV. 
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